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Abstract 


The reaction of HOC1 + HCI — > Ch + H 2 O in the presence of Cl has been studied using ab 
initio methods. This reaction has been shown to have a high activation barrier of 46.5 kcal mol' 1 . 
The chlorine anion, Cl is found to catalyze the reaction, viz. two mechanisms. The first involves 
Cl interacting through the concerted four-center transition state of the neutral reaction. The 
Other mechanism involves the formation of a HCTHOC1 • Cl intermediate which dissociates into 
CI 2 + Ci+ H 2 O. The steps are found to have no barriers. The overall exothermicity is 15.5 kcal 
mol' 1 . 


Understanding the underlying chemistry behind stratospheric ozone depletion is central to 
finding ways of curtailing further loss of ozone in the Antarctic Ozone Hole. One of the major 
reaction processes implicated in the catalytic removal of ozone is the heterogeneous reaction of 
HCI with CIONO 2 to give HNO 3 and CI 2. 1 ’ 2 ' 3 The CI 2 can subsequently photolyze to produce Cl 
atoms which can then participate in homogeneous catalytic cycles for ozone destruction through 

Cl + 0 3 -> CIO + 0 2 ( 1 ) 

CIO + O -4C1 + 0, ( 2 ) 

net: 0 + 03 —>202 (3) 

Recent field measurements of HCI and CIO concentrations in the stratosphere show that the 
reaction of HCI with CIONO 2 is critical in determining the winter and spring chlorine budget 
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within the polar vortex. 4,5 However, extensive experimental studies suggest that the 
heterogeneous reaction of HC1 with CIONO 2 , 

HC1 + CIONO 2 -4 HN0 3 + Cl 2 (4) 

may not occur in one step but two viz. ' ' 

ClONCh + H 2 0 -4 HOC1 + HN0 3 (5) 

HOC1 + HC1 + Cl 2 + H 2 0 (6) 

The net reaction of steps 5 and 6 is reaction 3. 

A key step is the reaction involving HOC1 with HC1 (reaction 6). Experimentally, 9 this 
reaction is known to be slow in the gas-phase; but little is known about the details of the 
energetic that could explain the slow rate. Abbatt et al. 6,7 have asserted that reaction 6 can occur 
on polar stratospheric cloud (PSC) particles. However, the existence of solvated chlorine in the 
form of Cion PSC particles has been suggested. 5,6,7 This has raised the question of whether Cl' 
can catalyze the reaction of HOC1 with HC1 (reaction 6). In this paper we use ab initio 
calculations to ascertain if Cl' can lower the activation barrier for the PIOC1 -t HC1 reaction. 

The geometries of reactants, products and transition states have been optimized using the ab 
initio MP2 method. 10 Two basis sets have been used for optimization: 6-3 lG(d) and 
6-31 1++G(2d,2p) which includes diffuses- and p-functions and an extra set of d-functions on all 
heavy atoms. The diffuse functions are added because we believe these to be essential for 
accurate description of geometries and energetic for anions. Single point energies have been 
calculated using the CCSD(T) method with the extend 6-31 1++G(3df,3pd) basis set. 

Geometries of the transition state for the reaction of HOC1 + HC1 —> Cl 2 + H 2 0 is shown in 
Figure la. The number without the asterisk are geometric parameters that have been calculated at 
the MP2/6-3 lG(d) level, and the numbers with the asterisk are at the MP2/6-31 1 ++G(2d,2p). 

The transition state is a four-center reaction. The IiCl approaches HOC1 orienting the hydrogen 


606 



of the HC1 toward the oxygen of HOC1. The HC1 bond length increases from the length of the 
HC1 bond in isolated hydrogen chloride. The CIO bond length also increases from that in HOC1. 
The transition state is a non-planar structure. The four atoms involved in the four-center ring is 
out-of-plane by 5.9°. Moreover, the hydrogen on HOClis out of the plane of the ClOH’ group by 
98°. The transition state is a true first-order saddle point, with one negative frequency. 11 The 
relative energetic are given in Table 1. The experimental heat of reaction 12 at OK for the HOC1 + 
HC1 reaction is -18.0 kcal mol' 1 . At the CCSD(T)/6-31 1++G(3df,3pd) level we predict a value of 
-15.4 kcal mol’ 1 , which is in reasonable agreement with the experimental results. The calculated 
barrier of the same level of theory is 46.5 kcal mol' 1 . Such a high barrier suggest that the 
homogeneous reaction will be very slow which is consistent with the experimental observations. 

Ln the presence of C1-, the anion prefers completing through the hydrogen not involved in the 
four-center ring. The anion has a dramatic effect on the transition state. The transition state 
along with structural parameters are given in Figure lb. The transition state is also a first-order 
saddle point with one imaginary frequency. 13 Comparing the HC1 bond, this bond is longer in the 
transition state involving the anion than without the anion present. The CIO bond in the transition 
is significantly elongated as a result of the C1-. It is also interesting to note that the atoms in the 
four center ring is more out-of-plane than in the neutral reaction; 24.0° versus 5.9°, respectively. 
These changes in the transition state for HOC1 + HC1 involving Cl compared to without the CT 
indicate that the transition state occurs earlier in the entrance channel. The structural changes in 
the transition state also suggest that breaking the CIO bond in HOC1 and the HC1 bond in HC1 in 
the transition state is easier by the presence of the C1-. The activation barrier for HOC1 + HC1 + 
Cl” formation is 8.3 kcal mol' 1 at the CCSD1TV6-31 1 ++G(3df,3pd) level. The activation barrier 
for the neutral reaction (3) is found to be very high and the reaction is slow. However, 
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involvement of Cl anion dramatically changes the potential energy surface. These calculations 
shows that the involvement of the anion lowers the neutral barrier by 38.2 kcal mol' 1 . Moreover, 
these calculations suggest that the anion could play a significant role in reducing the energetic of 
the HOC1 + HC1 reaction. 

While exploring the potential energy surface of the HOC1 + HC1 reaction in the presence of 
C1-, we found an alternate mechanism for the formation of the CI 2 + H 2 O products. A stepwise 
mechanism involving the formation of anion complexes such as Cl - HOC1 and Cl * H 2 O was 
found. The C1‘ ■ HOC1 anionic complex reacts with HC1 in a barrier-less process to form a HC1 - 
HOC1 . Cl intermediate. This intermediate is shown in Figure 2. At the MP2/6-31G(d) level, the 
CIO bond in the HC1 • HOC1 • Cl intermediate is elongated nearly 0.6 A from the CIO bond in 
HOC1, and the H'Cl bond is elongated by 0.7 A from the uncomplexed HC1. We have performed 
a vibrational frequency analysis to determine if the intermediate is a saddle or minimum energy 
point on the potential energy surface. 14 The vibrational energy calculations reveal that the 
structure has all real and positive frequencies. The energy of the intermediate formation from the 
reaction of HC1 + HOC1 and Cl is calculated to be -37.3 kcal mol' 1 at the MP2/6-3 lG(d) levels, 
and reduces to -38.0 kcal mol' 1 at the CCSD(T)6-311++G(3df,3pd) level, The formation of CI 2 + 
H 2 0 + Cl requires 22.5 kcal mol' 1 at the CCSD(T) level. As a result, one can view an alternate 
mechanism for HC1 + HOC1 catalysis by Cl' as involving a step-wise mechanism first involving 
the formation of a Cl . HOC1 complex which then attaches to HC1 forming an HC1 . HOC1 . Cl' 
intermediate. The next step involves the detachment of CI 2 and Cl . H 2 O complex. At low 
temperatures this mechanism involving the HC1 -HOC1 . Cl intermediate may be preferable. We 
should point-out that such an attachment/detachment mechanisms has been suggested for the 
reaction of HC1 with CIONO 2 catalyzed by N0 3 ‘. 


608 



More importantly, the present calculations suggest that anions, such as C1-, can play a major 
role in promoting the reaction of HOC1 with HCI. The energetic effects on the HOC1 + HC1 
reaction are considerable. Moreover, anions such as nitrate and/or sulfate which are important 
ions in atmospheric aerosols may play similar roles in the energetic of the HOC1 with HCI 
reaction. 
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U + HOC! + HU -> HU • HUU • U 

H0C1 + HC1 -> Cl 2 + H 2 0 [H0C1 + HC1 + Clf HC1 • H0C1 • Cl' Cl 2 + H 2 0 + Cl 
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Figure Captions 


Figure 1 . (a) Transition state for HC1 + HOC1 — > CI 2 + H 2 O for the neutral reaction. Numbers 
without asterisk are calculated at the MP2/6-3 lG(d) level, numbers with the asterisk 
are at the MP2/6-31 1++G(2d,2p) level of theory. 

(b) Transition state for HC1 + HOC1 + C1 + CI 2 + H 2 O + Cl-with Cl interacting 
through the four-center HC1 + HOC1 transition state. Numbers without asterisk are 
calculated at the MP2/6-3 lG(d) level, numbers with the asterisk are at the MP2/6- 
31 1++G(2d,2p) level of theory. 

Figure 2. HC1 • HOC1 • Cl intermediate. Numbers without asterisk are calculated at the MP2/6- 
3 lG(d) level, numbers with the asterisk are at the MP2/6-31 1++G(2d,2p) level of 
theory. 
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